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Abstract 

Recent second-order perturbation computations have provided an accurate prediction for the 
primordial gravitational potential, <&(x), in scenarios in which cosmological perturbations are gen- 
erated either during or after inflation. This enables us to make realistic predictions for a non- 
Gaussian part of <J?(x), which is generically written in momentum space as a double convolution 
of its Gaussian part with a suitable kernel, /NL(ki,k2). This kernel defines the amplitude and 
angular structure of the non-Gaussian signals and originates from the evolution of second-order 
perturbations after the generation of the curvature perturbation. We derive a generic formula for 
the CMB angular bispectrum with arbitrary /NL(ki, k2), and examine the detectability of the pri- 
mordial non-Gaussian signals from various scenarios such as single-field inflation, inhomogeneous 
reheating, and curvaton scenarios. Our results show that in the standard slow-roll inflation sce- 
nario the signal actually comes from the momentum-dependent part of /NL(ki,k2), and thus it 
is important to include the momentum dependence in the data analysis. In the other scenarios 
the primordial non-Gaussianity is comparable to or larger than these post-inflationary effects. We 
find that WMAP cannot detect non-Gaussian signals generated by these models. Numerical cal- 
culations for / > 500 are still computationally expensive, and we are not yet able to extend our 
calculations to Planck's angular resolution; however, there is an encouraging trend which shows 
that Planck may be able to detect these non-Gaussian signals. 

PACS numbers: 98.80.Cq, 95.35.+d, 4.62. +v 
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I. INTRODUCTION 



Inflation is a building-block of the standard model of modern cosmology. It is widely believed 
that there was an early stage in the history of the Universe - before the epoch of primordial 
nucleosynthesis - when the expansion rate of the Universe was accelerated. Such a period of 
cosmological inflation can be attained if the energy density of the Universe is dominated by 
;he vacuum energy density associated with the potential of a scalar field, called the inflaton 
l|. Inflation has become so popular also because of another compelling feature. It provides 
a causal mechanism for the production of the first density perturbations in the early Uni- 
verse which are the seeds for the Large-Scale Structure (LSS) of the Universe and for the 
Cosmic Microwave Background (CMB) temperature and polarization anisotropies that we 
observe today. In the inflationary picture, primordial density and gravity-wave fluctuations 
were created from quantum fluctuations and then left the horizon during an early period of 
superluminal expansion of the Universe, with the amplitude "frozen-in". Perturbations at 
the surface of last scattering are observable as temperature and polarization anisotropies in 
the CMB. The inflationary paradigm has been tested carefully by the data of the Wilkinson 
Microwave Anisotropy Probe (WMAP) mission [2j]. The WMAP collaboration has pro- 
duced a full-sky map of the angular variations of the CMB with unprecedented accuracy. 
The WMAP data confirm the inflationary mechanism as responsible for the generation of 
curvature (adiabatic) superhorizon fluctuations Q|. Since the primordial cosmological per- 
turbations are tiny, the generation and evolution of fluctuations during inflation has been 
studied within linear perturbation theory. Within this approach, the primordial density 
perturbation field is a Gaussian random field; in other words, its Fourier components are 
uncorrelated and have random phases. Despite the simplicity of the inflationary paradigm, 
the mechanism by which cosmological adiabatic perturbations are generated is not yet fully 
established. In the standard slow-roll scenario associated with one-single field models of in- 
flation, the observed density perturbations are due to fluctuations of the inflaton field itself 
when it slowly rolls down along its potential. When inflation ends, the inflaton oscillates 
about the minimum of its potential and decays, thereby reheating the Universe. As a result 
of the fluctuations each region of the Universe goes through the same history but at slightly 
different times. The final temperature anisotropies are caused by inflation lasting for differ- 
ent amounts of time in different regions of the Universe leading to adiabatic perturbations 
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Q. 

An alternative to the standard scenario is represented by the curvaton mechanism 
"s| where the final curvature perturbations are produced from an initial isocurvature 
perturbation associated with the quantum fluctuations of a light scalar field (other than the 
inflaton), the curvaton, whose energy density is negligible during inflation. The curvaton 
isocurvature perturbations are transformed into adiabatic ones when the curvaton decays 
into radiation, much after the end of inflation. 

Recently, other mechanisms for the generation of cosmological perturbations have been 



proposed, see H for a comprehensive review For instance the i„hon_ reheating see- 
nario [10|] acts during the reheating stage after inflation if superhorizon spatial fluctuations 
in the decay rate of the inflaton field are induced during inflation, causing adiabatic pertur- 
bations in the final reheating temperature in different regions of the Universe. Alternatively, 
curvature perturbations may be created because of the presence of broken symmetries during 
inflation 

Testing the Gaussianity of the primordial fluctuations provides a powerful tool to dis- 
criminate between different scenarios for the generation of the cosmological perturbations 
which would be indistinguishable otherwise . Non-Gaussianity is a deviation from a pure 
Gaussian statistics, i.e., the presence of higher-order connected correlation functions of CMB 
anisotropies. The angular n-point correlation function is a simple statistic characterizing a 
clustering pattern of fluctuations on the CMB. If the fluctuations are Gaussian, then the 
two-point correlation function specifies all the statistical properties of higher-order correla- 
tion functions, for the two-point correlation function is the only parameter in a Gaussian 
distribution. If it is not Gaussian, then we need higher-order correlation functions to deter- 
mine the statistical properties. For instance, a non-vanishing three-point function of scalar 
perturbations, or its Fourier transform, the bispectrum, is an indicator of non-Gaussian 
features in the cosmological perturbations. The importance of the bispectrum comes from 
the fact that it represents the lowest order statistics able to distinguish non-Gaussian from 
Gaussian perturbations. An accurate calculation of the primordial bispectrum of cosmo- 
logical perturbations has become an extremely important issue, as a number of present 
and future experiments, such as WMAP and Planck, will allow us to constrain or detect 
non- Gaussianity of CMB anisotropy with high precision. 

In order to compute and keep track of the non- Gaussianity of the cosmological perturba- 
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tions throughout the different stages of the evolution of the Universe, one has to perform a 
perturbation around the homogeneous background up to second order. Recent studies have 
been able to characterize the level of non-Gaussianitv predicted in the various scenarios for 
the generation of the cosmological perturbations 
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On large scales the second-order, gauge-invariant expression for the temperature 



anisotropies reads |9 
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where <3>l represents the gauge-invariant gravitational potential at linear order, is the lin- 
ear gauge-invariant comoving curvature perturbation, (2 is the second-order gauge-invariant 
comoving curvature perturbation, and 
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It shows that there are two contributions to the final nonlinearity in the large-scale temper- 
ature anisotropies. The third term, ((2 — 2Cl) 5 comes from the "primordial" conditions set 
during or after inflation. They are encoded in the curvature perturbation ( which remains 
constant once it has been generated. The remaining part of Eq. (0) describes the post- 
inflation processing of the primordial non-Gaussian signal due to the nonlinear gravitational 
dynamics, including also second-order corrections at last scattering to the Sachs- Wolfe ef- 
fect. Thus, the expression in Eq. (JTJ) allows us to separate the primordial contribution to 
non-Gaussianity from that arising after inflation. 

While the nonlinear evolution after inflation is the same in each scenario, the primor- 
dial term will depend on the particular mechanism generating the perturbations. We may 
parametrize the primordial non-Gaussianity in the terms of the conserved curvature pertur- 
bation (in the radiation or matter dominated epochs) (2 = 2onl (Cl) 2 , where onl depends on 
the physics of a given scenario. Within the standard scenario where cosmological perturba- 



tions are due to the inflaton the primordial contribution to the non-Gaussianity is given by 
«nl = 1 — \( n c ~ 1) 3 13, 3|, where the spectral index is expressed in terms of the usual 
slow-roll parameters as — 1 = — 6e + 2r/ [jj. In the curvaton case onl = (3/4r)— r/2, where 



r ps (prj P)~d is the relative curvaton contribution to the total energy density at curvaton 
decay |9j. In the minimal picture for the inhomogeneous reheating scenario, onl = 1/4. 

^From Eq. ((TJ) one can extract the nonlinearity parameter / NL which is usually adopted 
to phenomenologically parametrize the non-Gaussianity level of cosmological perturbations 



and has become the standard quantity to be observationally constrained by CMB experi- 
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2l\ . The comparison between our expression [Eq. (0)] and that in the previous 
2l| can be made through the Sachs- Wolfe formula, AT/T = —(1/3)$, where $ is 
Bardeen's gauge-invariant potential, which is conventionally expanded as (up to a constant 
offset, which only affects the temperature monopole) 

$ = $l + /nl*$l- (3) 

Here the ^-product (convolution) makes explicit the fact that the nonhnearity parameter 
has a non-trivial scale dependence Therefore, using ( L = — (5/3) $l during matter 
domination, from Eq. <$J$ we may define the nonhnearity parameter in momentum space 



/ N L(ki,k 2 ) 

where 



jj(l-a NL ) + ^-^/C 



(4) 



= ¥ 3F — ' (5) 



10 (ki • k 3 )(k 2 -k 3 ) _ 10_(ki ■ k 2 

and k 3 + ki + ko = and k 3 = | k 3 [ 1 . Notice that in the "squeezed" limit first discussed 
by Maldacena [13j, where one of the wavenumbers is much smaller than the other two, e.g. 
k\ /c 2 , 3 , the momentum dependence of the kernel disappears. 

The fact that the nonhnearity parameter has a scale (momentum) dependence, that is 
that /nl is not simply a number, may call for a reanalysis of the tests performed so far of 
the non-Gaussianity in the primordial cosmological perturbations I20I u\. This is because 
previous studies have been done when theoretical predictions for the nonhnearity parameters 
in the various scenarios (including the standard case in which perturbations are generated 
by the inflaton field) were not available and therefore /nl was assumed phenomenologically 
to be a constant. 

The observational capability of determining the nonhnearity parameter /nl is the sub- 
ject of a long project of which this paper represents the first step. Starting from a generic 
expression for the gravitational potential, we first derive the generic expression for the pri- 
mary CMB angular bispectrum. This formula generalizes the one provided by Komatsu 
and Spergel j20j| who worked with a constant /nl in momentum space. We then estimate 



1 The formula already accounts for an additional nonlinear effect entering in the CMB angular 3-point 
function from the angular averaging performed with a perturbed line-element implying a +1 shift in /nl- 
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the expected signal-to-noise ratio for detecting primary non-Gaussianity at WMAP angular 
resolution. While we show that the primary non-Gaussian signal generated in standard sce- 
narios of inflation cannot be detected by WMAP, our predicted signal-to-noise ratio shows 
a trend which, if maintained at higher angular resolution, should allow detection of the 
non- Gaussian signals by the future Planck mission even in the standard single-field scenario 
of inflation - in this case, /nl is dominated by the post-inflationary evolution, rather than 
the primordial contribution from inflation. 

The paper is organized as follows. In section |H] we give some basic definitions and we 
compute analytically the CMB angular bispectrum arising from a primordial potential of 
the kind described by equation P|): in section 11111 we present our numerical predictions 
for the primary angular bispectrum, and discuss detectability with the current and future 
experiments; section IIVI contains our concluding remarks. 



II. THE CMB ANGULAR BISPECTRUM 
A. Basics 

The CMB angular bispectrum is defined by 

^hMs 3 = { a t-\mi a i2m 2 a hmz) i (6) 

where we have expanded the observed CMB temperature fluctuations into spherical har- 
monics, and we have defined the multipoles ag m 

a im = J d n — - — Y lm . (7) 

We find it convenient to split the multipoles d£ m into a Gaussian part a\ m and a non-Gaussian 
part <^£ m ■ 

By ignoring second-order terms in a^, we obtain 




7 



The rotational invariance of the CMB sky implies that B™}^™ 3 can always be decomposed 

as 



Bhhh ■ (10) 



' ii i 2 k 
, nil m 2 m 3 

Where B^ 2 ^ 3 is the angle-averaged bispectrum and the matrix is the Wigner 3j symbol. 
The presence of the Wigner 3j symbol ensures that the bispectrum satisfies the selection 
rules, nii+m 2 + m 3 = 0, £i+£ 2 +£ 3 = even, and the triangle conditions, \£i— £j\ < £k < £i+£j 
for all permutation of indices i, j, k. 

As we have mentioned in the Introduction, in the various scenarios for the generation 
of the cosmological perturbations, the non-Gaussian part of the primordial gravitational 
potential can be expressed in Fourier space as a double convolution, 



$ NL (k 3 ) = J d 3 h Sk 2 5 (3) (kx + k 2 - k 3 ) $ L (ki) $ L (k 2 )/ NL (k 1 , k 2 , k 3 ) , (11) 

where $L(k) is a Gaussian random field representing the Gaussian part of the primordial 
potential; the kernel, /NL,(ki, k 2 , k 3 ), in equation (fTTj) can be written, without loss of gen- 
erality, as: 

/NL(ki, k 2 , k 3) = £ ^k^ki.hr . (12) 

n=0 fc 3 

in the following we are going to expand /NiXki, k 2 , k 3 ) in Legendre polynomials in terms of 
the angle between ki and k 2 2 : 

N 

/ NL (k 1; k 2 , k 3 ) = J2 Mb, h, k 3 )P e (k ■ k 2 ) . (13) 

The multipoles of the harmonic expansion of the (today observed) CMB temperature 
anisotropies are related to the primordial potential $(k), the relation between the two 
quantities being described by the linear radiation transfer functions, A^(&,t ): 



2 In equation we defined the kernel as a function of ki and k2 only, as k3 was given by k3 = ki + k2; 
nevertheless for our following derivation of the bispectrum we find it convenient to introduce the Dirac 
delta function 8^ (ki + k2 + k3) in eqn. and to write the convolution kernel as a function of kj, 

k2, k3 separately. In this way we can avoid factors of (k\ ■ A^)" in the denominator of i|12|) . Those factors 
would make the decomposition of the kernel in Legendre polynomials difficult. 
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aim = {-if J ^3 *(k) A e (k, r ) Y; m {k) , (14) 

where we are evolving the primordial perturbations up to the present time To- In the following 
we write simply Ae(k) instead of Ae(k, tq). 

The primordial potential is the sum of a linear and a nonlinear part: $(k) = $L(k) + 
$nl (k) , where the non-Gaussian part is given by formula (jllj) ; accordingly, we can split 
also the temperature fluctuation and the multipoles ae m into Gaussian and non-Gaussian 
components. Our aim in the next section will be to calculate the CMB angular bispectrum, 
starting from the bispectrum of the primordial gravitational potential which is, by definition 
($L(ki)$L(k2) <: &NL(k3)) + cyclic permutations. 

B. Analytic formula of the primary bispectrum with arbitrary kernel 



Let us first fix the notation by explicitly writing Eq. (fTlj) for a\ imv a\ m2 and a\ 

WHY'J^M^YUi^A^h) , (15) 



a £ imi 



a hm 2 



(47F)H) ' 2 / 1^ $L(k2) Y ^ ik) A ' M ■ (16) 

< L m3 = (4vr)(-i) 4 / 1^ $ NL (k 3 ) Y; 3m Xh) A £3 (h) . (17) 
Now, putting together Eqs. ([To]) . (fTBj) and (|T7j) . and using ©, we find 

*ssr = ( 4 -) 3 (- i )" +&+ '7Him < * L(ki) * i ' (k2) * Ni(k3)> 



. K 3 <-> ti 

+ 



(18) 



( 

The component ($L(ki) < I ) L(k2)$NL(k3)) of the <3>(k)-field bispectrum can be easily calcu 
lated: 

($ L (k 1 )$ L (k 2 )<I»NL(k3)) = 2(2vr) 3 ^ (kx + k 2 + k 3 )/ N L(k 1 , k 2 , k 3 )P(A; 1 )P(A; 2 ), (19) 
and we obtain 

Bissr = (^) 3 (-i)'' +&+ " / |^H02(2,)^<" (kl + k 2 + k 3 )/ NL( k 1 , k 2 , k 3) 
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x P(/ Cl )P(A; 2 )A 4 (A; 1 ) A h (k 2 ) A £a (k 3 ) 3£ mi (*i) Y e * 2m2 (k 2 ) Y^h) 
( h <-> £1 \ 1 1, ^£i\ 

+ 



( h <- h ^ 

























(20) 



The Dirac delta function 5^ (ki + k2 + k 3 ) can now be written as 

d 3 r 



6® (k x + k 2 + k 3 ) = J 



(2vr)= 



iki-r ik 2 -r ik 3 -r 



and the plane waves can be expanded according to the Rayleigh formula 

e* k ' x = (47r)J2J2{iYj e (kx)Y em (k)Y; m 



x 



(21) 



(22) 



l rn 



In this way we can make the substitution 



<S®> (kx + k 2 + k 3 ) = 8 £ J] S^" 4 *W^) y W^) y 4™' 3 (^) 

(23) 



where we have introduced the Gaunt integral G^i™?™ 3 , defined by 

dnY^ {n)Y^ m > 2 (n)Y e > 3m > 3 (n) 



y pi fi pi 



i' x £' 2 £' 3 



(24) 



, (2f 1 + i)(2f 2 + i)(2f 3 + i) ( e x e 2 e 3 





An 



( p> pi pi 

l l l 2 c 3 

\ m' 1 m' 2 rn' 3 



The kernel /NL,(ki, k 2 , k 3 ) can be expanded in spherical harmonics as well: Eq. (|13|). together 
with the addition theorem of spherical harmonics, 

4tt l 



P L (h-k 2 ) = — [ £ Y LM (h)Y£ M (k 2 ) 



(25) 



M=-L 



finally yields 



N 



/ NL (ki,k 2 ,k 3 ) = h{h, k 2 ,k 3 )P L (ki ■ k 2 ) 

L=0 



N 



Air 



E^—jfL(knk 2 ,k 3 ) £ Y LM (ki)Y£ M (k 2 ) 



L=0 M=-L 

Now, splitting the integral on the right-hand side of Eq. (J20|) into a radial and an angular 
part, and considering Eqs. (J25j) and (J25J) . we find 



7T 



( 2L + !) 



,2 rLt'A 
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E 

Mm!^vn! 2 



-i \mi+mi pm\m' 2 m:j „-m 1 Mm' l ( -,-m' 2 m2M 

~ L ) y e , 1 &' 2 £ 3 ^iii; ^£ 2 i 

































« 3 J 



+ 



where we have used orthonormality of spherical harmonics, and we have defined 



x 



dk 2 klP< s ,(k2)A e2 (k2)ji>Ak 2 r)f L (k 1 , fa, k 3 



(27) 



(28) 



Formula ()27|) is what we have been looking for: it describes the angular CMB bispectrum 
arising from the primordial potential [Eq. (fTTj) ]. The angle-averaged bispectrum, B^ 2 ^ 3 , is 
related to B™}^™ 3 by Eq. (fTUj). and an explicit expression for the angle-averaged bispectrum 
can be easily derived from Eq. ()27|). We use the following relation of the Wigner symbols: 



E (- 

Mm^ m' 2 



mi +m' 2 pm' x m' 2 m 3 ~—m\Mm' x ^-,-m' 2 m2M 

y i' x t' 2 e 3 ^i l e[ y i' 2 1 2 l 



^ «2 «3 
? 2 ^'i 



X 




(29) 



'1 (2 13 

where < ^ is the Wigner Qj symbol, and we have defined the quantities 

? 2 £ l L 



l iie 2 is 



'(24 + 1) (2£ 2 + l) (2£ 3 + l) ( t\ h 





47T 



(30) 



Using these quantities, we obtain the final analytic formula of the angle-averaged bispectrum 
with arbitrary kernels: 



B 



■ ' L= i>i> 2=Q 2L + 1 I £' 2 £[ L 



EE 



x drr 2 C f JJ 2 (r) + 



-3 w t-2 
?1 4 



(31) 



'2 w t-3 

We use this general relation to calculate numerically the CMB angle-averaged bispectrum 
for the class of inflationary models that produce potentials in the form of Eq. (J11J) . To select 
and study a specific model we need to provide an explicit expression for the coefficients of the 
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Legendre expansion of the kernel, Eq. (|13p (i.e. we need to provide an explicit expression 
for &2, ks) in Eq. (|2*%j) ). We will now consider the various possibilities for the kernels 

in the next subsections. 



C. Constant Kernel 



The simplest possible choice of the kernel is a constant, /NL,(ki, k 2 , k 3 ) = / NL (where /nl is 
a constant parameterizing the level of non-Gaussianity), which gives 



<&(x) = $ L (x) + /, 



NL 



d>*(x) - (<^(x) 



(32) 



in real space. This is the usual phenomenological parametrization of non-Gaussianity which 
has been widely used in the literature. The CMB angular bispectrum in this model has been 
calculated by Komatsu and Spergel j^. As a simple check of our calculations, we rederive 
their formula starting from Eq. (|31jl . 

For a constant /nl, N = and fo(ki, &2, ks) = /nl in Eq. thus, Eq. (|3*T|) yields 



2 oo 



'1 *2 ^3 



~ 7 - ~ ' * £ o 



44=0 



We can write 



h' 2 t 2 o 



-3 w 
?2~4 

' (2*2 + 1) , 
4vr 1 



where 5 J is a Kronecker delta and we have used the formula: 



4 4 

m — m 



r 



An analogous relation holds for I^o, giving 

R /8\ 2 u e 1+h+t3 (2l 3 + l)§( 2l 2 + l)(2l 1 + l) 



(47T) 



X 



-1 t-2 t 3 
^2 4 

„2 r 0£^ 2 



Kl *2 *3 





^4^4^ 










^4 <- 4 1 








v 4~4 J 




^4^4 J 



(33) 



(34) 



(35) 



(36) 
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If we now make use of the relation 



'1 t-2 t 3 
? 2 £l 



-1) 



v /(2£ 2 + l)(2£ 1 + l) 
and remember that, when / NL (ki, k 2 , k 3 ) = / NL , by definition (see Eq. 



(37) 



/nl / dk 1 kl^i x {kx)P^{kx)jt- L {k 1 r) \ dk 2 klP0(k 2 )A h (k 2 )ji 2 (k 2 r) 



x / dk 3 klA h (k 3 )j h (k 3 r) 
we recover exactly the result of Ref. 



(38) 



D. Momentum-dependent kernels 

Throughout the rest of this work we are going to consider a primordial potential kernel 
defined by Eq. (gj): 



/ NL (ki,k 2 ,k 3 ) 



1 5 
6~3 



fci • k 2 , kfk% t ou2u2 (h ■ k 2 ) 
'A 

v 3 



(1 - Qnl) k\k 2 "\ 2 "' +3'- z jj £ +3k'ik' 2 : — 4 "" +3fcifc 2 

/Co /Co /Co 



(^ + fc 2 2 )( fcl . fc2 ) 



^ : 3 



(39) 

It follows from the form of the kernel that we can expand /NL,(ki, k 2 , k 3 ) into the first three 
Legendre polynomials in terms of the angle between k x and k 2 : 



/ NL (k l5 k 2 , k 3 ) = k 2 , k 3 )Pt(k! ■ k 2 ) 



PoCh ■ k 2 ) 

Pi(h-k 2 ) 
P2(k ■ k 2 ) 



Pi(fci • k 2 ) = h-k 2 , 



A simple, direct calculation shows, for our kernel, that 



fo{ki,k 2 ,k 3 ) 
fi(h,k 2 ,k 3 ) 
f 2 (h,k 2 ,k 3 ) 



-i-jj(l-flNL) 



+ 



Akjkj 
k 3 



k\k 2 
k 3 

^3 



(40) 

(41) 
(42) 
(43) 

(44) 
(45) 
(46) 
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Therefore, we find that the conventional momentum-independent parameterization, /nl, 
captures only the first term in /q. We evaluate numerically the expression of the CMB 
angle-averaged bispectrum, which is obtained by substituting these fe(ki, k 2 , k 3 ) coefficients 
into (J2HJ): 



£ 



0£[£' 2 

e 3 £i£2 1 



(1 - Q-NLy 



^J(rWA( r W 2 k(r) +Hrmt'e 1 (rW 2 k(r) , (47) 



?(0) 



?(0) 



}(2) 



5(2) 



,u[f 2 



-«t a) (r)^(r)^(r) + 3a^(r)^(r)^(r) 
+3«t 4) (r)/3S 1 (r)^,(r) , 



(48) 



4ttW = 2aJ 3 - 4) (r)^S 1 (r)^(r) ; 



the quantities a and (3 being defined as 



dk k 2 k ni AfXk)j e (kr) 



(49) 



(50) 



/3g(r) = Jdhklk^P^A^j^hr) 



(51) 



We then use these results in Eqs. (J3T|) to compute the angle-averaged bispectrum numerically. 



III. NUMERICAL RESULTS 



A. Radial Coefficients 



The problem of the numerical evaluation of B^£ 2 £ 3 can be divided into two parts. The first 
is the calculation of the Wigner 3j and 6j coefficients, while the other is the generation of 
the coefficients £^^(r). Since the expansion of our kernel contains only the first three 
Legendre polynomials, we consider only < L < 2. This allows us to use analytic formulae 
of the Qj symbols. Also the 3j symbols in h^ 2 t 3 can be evaluated by the well-known analytic 
formulae based on the Stirling approximation at high ts. 

The calculation of £^^(r) can be reduced to the numerical evaluation of a\ ni \r) 
(Eq. [HO] ) and f3i^ 2 2 \r) (Eq. in which we have to account for all the possible choices 
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of the set of values {L , £\ , £ 2 , £3 , £[ , £' 2 }, while for rij we need only n\ = —4, —2, and 0, 

and n 2 = 0, 1, 2, and 3. (See Eq. Eq. (}3*T]) . applied to our case, shows that if we 

want to calculate a particular mode of the averaged bispectrum, B^^, we have to gener- 
ic i' £' 

ate all the terms of C^l^yr) for < L < 2, 1 < £[, £ 2 < 00. The selection rules of the 
Wigner coefficients guarantee that the only terms which contribute to the sum over £[, £' 2 
(for fixed L) are those of which satisfy the triangular conditions: £\ — L < £[ < £\ + L and 
£ 2 -L<£' 2 <£ 2 + L. 

In our analysis, we consider a Concordance Model with fl\ = 0.7, Qb = 0.05, fl m = 0.3, 
h = 0.65 and n = 1; Figures (JTJ) and (j2J) show some radial coefficients oq (r) and /3^" (r), 
calculated at the time of decoupling, r* = c(r — r*), where r denotes conformal time, r is 
the present conformal time, and the decoupling time, r*, is defined at the peak of visibility 
function. In our model we have ctq = 14.9 Gpc and cr* = 289 Mpc. To calculate the radial 
coefficients we use a modified version of the CMBfast code 3 . 

Although most of the signal is generated in a narrow region around decoupling (i.e. 
when r ~ r* = c(tq — r*) ~ ctq), in the low-£ regime we still have to account for the low-r 
contribution due to the late integrated Sachs- Wolfe effect. Thus our r- integration boundary 
is c(r — 6t*) < r < cr for £ > 50, whereas < r < cr for £ < 50. The step-size 
Ar is determined by the ratio of the width of the last scattering surface to the present 
cosmic horizon, cr , and by the necessity of an accurate sampling of the acoustic oscillations 
at recombination. As the number of oscillations increase at high-£, we need smaller and 
smaller step sizes when simulating experiments with higher and higher angular resolutions. 

B. Signal-to- noise ratio for WMAP 

Even if a significant angular bispectrum was detected in CMB, this would not necessarily 
mean that it was generated by some primordial mechanism like inflation. There are in fact 
several foregrounds which can produce non-Gaussianity in CMB anisotropics, such as the 
Sunyaev-Zel'dovich (SZ) effect, weak lensing, the presence of point sources, and so on. Thus 
a complete study of detectability of the primary bispectrum needs to include the secondary 

In our numerical computation we are neglecting second order corrections to the CMB radiation transfer 
functions. These should in principle be included for a complete and definitive treatment of CMB non- 
Gaussianity. 
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FIG. 1: Radial coefficients a ™ (r) [Eq. (|50|)] at the time of decoupling, r*. From top to bottom we 
plot (£- !)£(£ + !)(£ + 2)o^ 4) (r,), + l)a[~ 2) (r*), aj 0) (r*), respectively. 



bispectra generated by the foregrounds, in order to check if the primordial component can 
be isolated from others. 

Having calculated numerically the angle-averaged bispectrum from primary and sec- 



ondary sources, we evaluate a x 2 statistic 20] 



2 \ - \Ptt2i3 A i B eXi 3 

x = L -2 > ( 52 ) 



2<e 1 <£ 2 <e 3 a hhh 



where B^f 2i is the observed bispectrum and B^^ are the theoretically calculated bispectra 
for different components, denoted by i. The variance v\ x ^ 3 of the bispectrum can be written 
as 22, 3 
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FIG. 2: Radial coefficients @\ {r) [Eq. (|51j)] at the time of decoupling, r*. On the left side, from 
top to bottom: + l/2)/3^(r*), + 1)/% (r*). On the right side, from top to bottom: /?if ( r *)s 



a 



2 



B hhh) ~ ( B £iM 3 ) — C^C^C^A^n^ , (53) 



where takes values 1, 2, and 6 when all €s are different, two of them are equal and all 

are the same, respectively. Ce is the sum of the theoretical CMB angular power-spectrum 
and the power-spectrum of the detector noise. The last one can be calculated analytically 



using Ref. 



24|. 



Taking dx 2 /dAi = 0, the Fisher matrix is given by 
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n(i) d(j) 

F .. = £ ^|^3 ; (54) 
and the signal-to-noise ratio (S/N)i for a component i is 

V M 

Let us neglect for the moment the non-diagonal components of the Fisher matrix; then, 
denoting the primordial component by % — 1, we can give an estimate of the expected signal- 
to-noise ratio for the primordial non-Gaussian signal without considering foregrounds. It is 
simply 

'Fn ■ (56) 



We have calculated the approximated signal-to-noise ratio using formula (JBTjj) for an ex- 
periment with the FWHM beam-size of WMAP (FWHM = 13', £ ma x = 500), assuming 
different scenarios for the generation of the cosmological perturbations, namely the stan- 
dard single-field slow roll scenario, the inhomogeneous reheating scenario and the curvaton 
scenario. The shape of the kernel in all these scenarios is given by Eq. ((HSj) with model- 
dependent values of the constant onl- According to Ref. jjj, in single-field slow roll inflation 
onl = 1, in the inhomogeneous reheating case onl = 1/4, whereas in the curvaton scenario 
we have onl = (3/4r) — r/2, where r is the relative curvaton contribution to the total energy 
density at curvaton decay. 

Let us now comment on our results, starting from the standard single-field inflation and 
the inhomogeneous reheating cases. Even though we ignore correlations between the primor- 
dial and secondary bispectra, we find, for the standard inflationary scenario, the expected 
signal-to-noise ratio for WMAP is S/N ~ 0.10 and, for the inhomogeneous reheating case, 
S/N ~ 0.15; thus, the primordial non-Gaussianity from these models is below the WMAP 
detection threshold. As the correlation between the primordial and secondary bispectra 
would only lower the signal-to- noise ratio for the primordial component, we can conclude 
that the primordial bispectrum from these scenarios is undetectable with WMAP. 

These expectations confirm the ones obtained in the previous work where the non- 
Gaussian primordial gravitational potential was approximated as $nl = ^l + Znl^l an d /nl 
is a momentum-independent parameter defining the level of predicted non-Gaussianity. In 
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this framework Komatsu and Spergel derived a detection threshold of /nl = 20 for WMAP 
and /nl = 5 for Planck. This suggested that a primordial signal from standard single-field 
inflation would be undetectable by WMAP as, in this phenomenological approach, / NL was 



expected to be ~ 1 in the standard scenario 
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26, 
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C. Challenges of numerical calculations at high I 

Before we study the expected signal-to-noise ratio for the future high-resolution exper- 
iments such as Planck, we note that the kind of computation we have described so far is 
numerically very challenging. Even after parallelizing and optimizing as much as possible, 
our algorithm {e.g. by implementing analytic approximations for the Wigner coefficients 
and by minimizing the number of points in the integration samples), the highest ^max we 
can reach is £max = 500, corresponding to the angular resolution of the WMAP satellite. 
We have not been able to go beyond the WMAP resolution, as the CPU time requirement 
was too demanding. The parallelized version of our algorithm took 6 hours on 100 proces- 
sors to calculate the full bispectrum up to ^max = 500. As the CPU time scales roughly 
as ^max ^ would take about 5 years on the same number of processors to calculate it up 
to ^max — 3000, thus making approximations necessary. On the other hand, extrapolating 
our results to higher angular resolutions suggests that the primordial non-Gaussian signal 
could be significant enough to allow detection of the primary signal at £ ma x = 3000 (the 
angular resolution of the Planck satellite), even in the most standard single- field slow roll 
inflationary scenario. This will be explained in more detail in the next section. 



D. Prospects for detecting non-Gaussianity by Planck 

Komatsu and Spergel pointed out that even an ideal experiment needed /nl > 3 
in order to detect primordial non-Gaussianity. This last statement, when combined with 
the previous theoretical expectations of the amplitude of non-Gaussianity, /nl — 1, implied 
rather pessimistic prospects for detecting the primordial non-Gaussian signals in standard 
scenarios of single-field inflation. 

However, we stress here that the previous expectation, /nl — 1, even though it roughly 
took into account the effect of the post- inflationary evolution of non-Gaussianity, was not 
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based on the detailed second-order computation of the cosmological perturbations during 
and after inflation. For this reason, it must be considered only as an order of magnitude 
estimate, and care must be taken when we study detectability of non-Gaussianity from 
standard single-field inflation by the future experiments at high angular resolution such as 
Planck. 

Our prediction based on the complete second-order calculation of the primordial gravi- 
tational potential shows an encouraging trend which shows that the actual signal-to-noise 
ratio is larger than the previous prediction with /nl = 1- Even though it is insufficient to 
push the primary signal over the detectability threshold of WMAP, it could be big enough 
to allow detection of the primordial non-Gaussianity signals by Planck. 

Let us elaborate on this point. We evaluate the signal-to-noise ratio for 10 < ! ma x < 500 
both in our cases and in the standard /nl parametrization. In the standard parametrization 
of non-Gaussianity the signal-to-noise ratio can be written as 4 : 



where -^^(/nl = 1) is the Fisher matrix of the standard /nl model with /nl = 1- The 
idea is as follows: by comparing the actual signal-to-noise ratio predicted from our full 
calculations, (S/N) iul1 , to the standard one, we can estimate /nl that is required to produce 
the same (S/N) sta,nd in the standard parametrization as (S/N) : 



which is the /nl that is needed in the usual parametrization of non-Gaussianity to reproduce 
the same level of non-Gaussianity predicted by our model for a given £ max - This parameter 
allows us to compare the previous estimates to ours more easily. The results are shown in 
figure El where we consider two experiments with the beam and the noise characteristics 
similar to WMAP and Planck. Two things are worth noticing. First of all, /^{/(^max) is not 
constant over £ max . Second of all, / N i/(^max = 500) is significantly bigger than the previously 
expected value, /nl — 1, though it is still of the same order of magnitude. When we look 
at the Planck experiment, we also notice that f^{f(£ m ax) is monotonically increasing when 
^max > 40, reaching a value of / N i/(^max = 500) ~ 4,which is already very close to the 

4 We ignore the contributions from the foregrounds. 




stand 



(57) 




(58) 
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FIG. 3: Values of f^{ f (£ max ) (Eq. [SHI) in the standard single-field and inhomogeneous reheat- 
ing scenarios. This parameter, /^l (^max)> represents /nl in the usual parametrization of non- 
Gaussianity to reproduce the same level of non-Gaussianity predicted by our model for a given 
f mM . For the standard single-field inflation the contribution to non-Gaussianity comes only from 
the post-inflation non-linear processing of perturbations, which is independent of the inflation- 
ary model. Thus the solid line in the plots also represents that part of the non-Gaussian signal 
which must be present in the CMB anisotropics, regardless of the considered inflationary model. 
The lower panel shows our results for an experiment with beam size and noise characteristics 
similar to WMAP. The upper panel shows the same analysis for Planck. We are considering 
^max = 500, corresponding to the angular resolution of WMAP. A full analysis for Planck would 
require f m ax = 3000, which is beyond the current computational power (see Sec llll Qj) . 
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detection threshold /nl = 5 computed by Komatsu and Spergel [20j for the full resolution 
of Planck. Considering that the angular resolution of the Planck satellite corresponds to 
^max = 3000, and we stopped our computation at £max = 500, our results suggest that the 
non-Gaussian signal from standard single-field inflation is likely to be detected by Planck. 

In addition to the standard single-field inflation and the inhomogeneous reheating models, 
we also investigate the curvaton scenario. In this case the value of onl depends on the 
parameter r, the relative curvaton contribution to the total energy density at curvaton decay, 
as previously pointed out. We consider different values of r, and calculate /^L^max)- Note 
that the momentum- independent part has been calculated as /nl = — 5/(4r) + 5r/6 jl^ . 
Our results are summarized in table We notice that, for small values of r, the parameter 
/nl (^max = 500) is now smaller than what was expected in the previous predictions (meaning 
that our signal-to-noise ratio is smaller than what was predicted assuming the standard /nl 
parametrization), whereas f^{f{i max = 500) > /nl for r > 0.5. Therefore, it is incorrect to 
conclude that the amplitude of non-Gaussianity is smaller for larger r; on the contrary, the 
signal-to-noise stays nearly the same regardless of r. 

Before concluding this section, let us stress again that our rough estimate of / N {/(4nax) 
and the extrapolation of our results to the angular resolution of the Planck satellite does not 
allow any conclusive statement about detectability of the primordial non-Gaussian signals 
generated by the simplest models of inflation. It is important to keep in mind that our 
/nl (Anax) depends on £ max ; thus, it is still possible that f^{f{i max ) might start to decrease 
for £ max > 500 and stay below the detection threshold of Planck at £ max = 3000. However, 
we find an encouraging trend that / N { (4nax) increases monotonically for ^ max > 40. It is 
certainly worth finding a way to achieve the full numerical computation of the primordial 
bispectrum (|3Tj) at very high multipoles (^max — 3000). The current algorithm based on 
the full numerical integration of equation (|3Tjl is computationally very expensive, and most 
likely some approximations must be invoked; one possibility is to implement the flat sky 
approximation at high ts. This will be the topic of a forthcoming publication j^j]. 

IV. CONCLUSIONS 

In this paper, we have shown that the full second-order calculations of cosmological per- 
turbations and inflationary dynamics suggest that the realistic form of non-Gaussianity, 
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r 




feff 

NL 


S/N 


0.1 


12.42 


8.42 


0.24 


0.2 


6.10 


2.98 


0.08 


0.3 


3.92 


2.62 


0.06 


0.4 


2.79 


2.58 


0.07 


0.5 


2.10 


2.98 


0.08 


0.6 


1.58 


3.31 


0.10 


0.7 


1.20 


3.60 


0.10 


0.8 


0.90 


3.84 


0.10 


0.9 


0.64 


4.04 


0.11 


1.0 


0.42 


4.22 


0.12 



TABLE I: Results for the curvaton model. The first two columns show the value of the relative 
curvaton contribution to the total energy density at curvaton decay and the predicted values of 
/nl in the previous parameterization which assumes that /nl is a constant. The last two columns 
contain the new computation of Z^r (4nax = 500) (Eq. |58| ) and of the signal-to- noise ratio for 
WMAP. 



the kernel / N L(ki, k 2 , k 3 ), must contain momentum-dependent terms. We have derived the 
analytic formula for the angle-averaged primary CMB angular bispectrum. This formula 
allows a more realistic description of non-Gaussian CMB anisotropy, extending the phe- 



nomenological model adopted in Ref. 20], where /nl was taken to be a constant. We have 
developed a numerical code to compute the primary bispectrum and estimated the expected 
signal-to-noise ratio for detecting primary non-Gaussianity at the WMAP angular resolu- 
tion. Our results show that, in the framework of standard single-field inflation, the primary 
non-Gaussian signal cannot be detected by WMAP, as already indicated by the previous 
analysis. On the other hand, in our complete second-order approach to perturbations during 
and after inflation, we have found that the previous theoretical expectation, /nl — 1, was 
too pessimistic, and the actual value which defines the CMB bispectrum is much larger. This 
result implies that the primordial non-Gaussian signals might be detectable by the future 
Planck mission even in the standard single-field scenarios of inflation. However, using the 
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current numerical algorithm, we have not been able to reach Planck's angular resolution, 
Cax = 3000, which would require 5 years of CPU time on 100 processors, and our conclusion 
on the prospect for detecting non-Gaussianity by Planck has to rely on extrapolations from 
4ax = 500. Suitable approximations at high £'s will be required in the future, in order to 
make a definitive conclusion on detectability of primordial non-Gaussianity in CMB. Finally, 
let us comment on statistical methods to measure the bispectrum. Komatsu, Spergel and 
Wandelt have shown that the direct measurement of all possible configurations of the 
bispectrum is computationally too expensive, and developed a faster estimator of /nl as- 
suming that /nl is a constant. Recently, Creminelli et al. jsjj] have extended this method to 
the case where the dominant signals come from the equilateral configurations, which yields 
a certain momentum-dependence in /nl- Their model (Eq. [14] of |30j]), however, is different 
from the form of /Wi^ki, k 2 , k 3 ) in equation (|3U|) . and thus their estimator cannot be used 
to measure primordial non-Gaussianity from second-order perturbations. New estimators 
optimized to our /Ni^ki, k 2 , k 3 ) need to be developed. 
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